Comment on "Ising model on a small world network 
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In the recent study of the Ising model on a small-world network by A. P§kalski [Phys. Rev. E 64, 
057104 (2001)], a surprisingly small value of the critical exponent (3 w 0.0001 has been obtained for 
the temperature dependence of the magnetization. We perform extensive Monte Carlo simulations 
of the same model and conclude, via the standard finite-size scaling of various quantities, that the 
phase transition in the model is of the mean-field nature, in contrast to the work by A. Pekalski but 
in accord with other existing studies. 
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A. Ptjkalski in Ref. [jD studied the Ising model on 
a small- world network constructed from a ring lattice. 
In the presence of a finite fraction of additional long- 
range interactions, the model was observed to undergo 
a phase transition at a finite temperature, in agreement 
with other related studies ^ |[ Q. However, the ex- 
ponent (3, describing the critical behavior of the magne- 
tization in the vicinity of the transition, was found to 
be very small: (3 ~ 0.0001, in contrast to the previous 
studies suggesting the mean-field nature of the transi- 
tion (5[ ||, [|. The small- world network in Ref. Q was 
constructed in a slightly different way compared with the 
original model by Watts and Strogatz [|| , under the ad- 
ditional constraint that not more than one shortcut is 
allowed for each vertex in the network. In this comment 
we present results of extensive Monte Carlo (MC) simu- 
lations of the same model as that in Ref. [|, which reveal 
that the phase transition is described by the mean-field 
exponents, a — 0, (3 — 1/2, 7 = 1, and v — 2 (see below 
for definitions), and thus conclude that the additional 
constraint in the network construction does not change 
the universality class of the transition. 

For simplicity, we consider only the network where ev- 
ery vertex has one shortcut (or every spin has three cou- 
plings), which is identical to Pekalski's network with the 
parameter p = 1 in Ref. [!} We then perform extensive 
Monte Carlo simulations of the Ising model described by 
the Hamiltonian: 
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where J is the coupling strength, <7j(= ±1) is the Ising 
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spin on vertex i, and Aj denotes the neighborhood of ver- 
tex i, including those vertices connected to i. For given 
network size at temperature T (in units of J/Ub), we 
have measured various thermodynamic quantities such as 
Binder's cumulant 0, the specific heat, and the suscep- 
tibility: 
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with m = \(1/N) where (■ • ■ ) and [•■•] repre- 

sent the thermal average (taken over 5000 MC steps af- 
ter discarding 5000 MC steps for equilibration at each 
temperature) and the average over different network re- 
alizations (taken over 400-1200 different configurations), 
respectively. 

Binder's cumulant in Eq. (Q), plotted for various sizes, 
yields a unique crossing point as a function of the tem- 
perature T, providing a convenient method to determine 
the critical temperature T c . Figure |l](a) shows that the 
result T c = 1.81(2) is obtained from the crossing point 
of Binder's cumulant for A^ > 800. The critical exponent 
v, describing the divergence of the correlation volume in 
such a way that £y ~ \T — T c \~ u can be determined 
from the expansion of Um near T c 
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where T\ and T 2 (> T\) are chosen near T c . Figure |1|( 
results in the value v i=s 2.0, which, together with the hy- 
perscaling relation v = 2 — a, gives the critical exponent 
a w for the specific heat C v . 

With such a mean-field value, we write the finite-size 
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FIG. 1: (a) Binder's cumulant Un has a unique crossing point 
at T c — 1.81(2) (in units of J/fca). (b) The critical exponent 
v = 1.99 is obtained from the least-square fit to the form 
in Eq. (^). (c) Specific heat C v (in units of ka) also has a 
crossing point at T c — 1.82(2), suggesting a — 0. (d) From 
the similar expansion of C v near T c , v = 2.02 is obtained, 
(e) Finite-size scaling of the susceptibility again determines 
Tc = 1.82(2) with the critical exponent 7 = 1, (f) Finite- 
size scaling of the magnetization [(m)] [see Eq. (^)] with /3 = 
1/2 also leads to the crossing point at T c = 1.82(1). In (a), 
(c), (e), and (f), simulations have been performed with the 
temperature increment AT = 0.005 whereas the data for N > 
200 have been used for fitting in (b) and (d). 



scaling in the form 

C v = f^T-T c )N 1 ^. 



where f(x) is an appropriate scaling function with the 
scaling variable x. As shown in Fig. |l](c), the unique 
crossing point of C v again yields T c — 1.82(2), in ac- 
cord with T c obtained from Un within numerical errors. 
The similar expansion of C v then provides an alternative 
way of determining D: AC V = C V {T X ) - C V (T 2 ) oc N 1 /", 
leading to the estimation v w 2.0 in Fig. 0(d). 

The divergence of the susceptibility when T c is ap- 
proached from above is described by the critical exponent 
7: x ~ {T ~ r c )~ 7 , which suggests the finite-size scaling 
form 
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with the appropriate scaling function g(x). Combined 
with v ?» 2 found above, the finite-size scaling form (Q) 
yields the value 7 w 1 and T c = 1.82(2) as shown in 
Fig. 0(e). 

Finally, on the basis of the above observation, the crit- 
ical exponent (3 for the magnetization is then determined 
from the scaling form 



[(m)] = N-^ 9 h((T-T c )N 
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which leads to T c = 1.82(1) and f3 w 1/2 [see Fig. |(f)]. 

We note that the obtained critical temperature aj: 
pears to be higher by factor two than that in Ref. 
The standard mean-field approximation applied to this 
model, where the coordination number is three, yields 
Tmf = 3 1. Considering the infinite range of the in- 
teractions and the mean-field nature of the transition, 
we believe that our estimation T c = 1.82(2) (in units of 
J/fcfl) is more precise. In this respect, it is interesting 
to note that also in the XY model on the small-world 
network a relatively large value T c /Tmf ~ 0.9 has been 
estimated Q. 

In summary, we have numerically studied the Ising 
model on the small-world network constructed in the 
identical way as in Ref. 0. In contrast to Ref. [l| we 
have obtained the standard mean-field critical exponents: 
a = 0, p = 1/2, 7=1, and v — 2 and confirmed that the 
phase transition is of the mean-field nature, in agreement 
with other previous studies (2[ [| [| . 
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